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j-H , Abstract 

' Recent observations on ultra high energy cosmic rays (those cosmic rays with energies greater 

than ~ 4 X 10^** eV) suggest an abundant flux of incoming particles with energies above 1 x 10^° eV. 
These observations violate the Greisen-Zatsepin-Kuzmin cutoff. To explain this anomaly we argue 
that quantum-gravitational effects may be playing a decisive role in the propagation of ultra high 
energy cosmic rays. We consider the loop quantum gravity approach and provide useful techniques to 
establish and analyze constraints on the loop quantum gravity parameters arising from observational 
data. In particular, we study the effects on the predicted spectrum for ultra high energy cosmic rays 
and conclude that is possible to reconcile observations. 
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O ■ 1 Introduction 

. In the present discussion we are concerned with the observation of ultra high energy cosmic rays (UHECR), 

' i.e. those cosmic rays with energies greater than ~ 4 x 10^* eV. Although not completely clear, it has 

been suggested that these high energy particles are possibly heavy nuclei [1] (we will assume here that 
^ ^ they are protons) and, by virtue of the isotropic distribution with which they arrive to us, that they 

O . originate in extragalactic sources. 

The propagation of UHECR in intergalactic space is subject to interaction with the cosmic microwave 
background radiation (CMBR). Its presence produces friction on UHECR making them release energy in 
the form of secondary particles. This affects their possibility to reach great distances. A first estimation of 
the characteristic distance that UHECR can reach before losing most of their energy was simultaneously 
' made in 1966 by K. Greisen [2] and G. T. Zatsepin & V. A. Kuzmin [3]. They showed that the observation 

of cosmic rays with energies greater than 4 x 10"'^^ eV should be greatly suppressed. This energy is usually 
referred as to the GZK cutoff energy. A few years later, F.W. Stecker [4] computed the mean life time for 
protons as a function of their energy, giving a more accurate perspective of the energy dependence of the 
cutoff and showing that cosmic rays with energies above 1 x 10^° eV should not travel more than ~ 100 
Mpc. More detailed approaches to the GZK-cutoff feature have been made since these first estimations. 
For example V. Berezinsky & S.I. Grigorieva [5], V. Berezinsky et al. [6] and S.T. Scully & F.W. Stecker 
[7] have made progress in the theoretical study of the spectrum J{E) (i.e. the flux of arriving particles 
as a function of the observed energy E) that UHECR should present. As a result, the GZK cutoff exists 
in the form of a suppression in the predicted flux of cosmic rays with energies above ~ 8 x 10^^ eV. 

Currently there are two main different sets of data for the observed flux J{E) in its most energetic 
sector [E > A X 10^^ eV). On one hand we have the observations from the High Resolution Fly's Eye 
(HiRes) collaboration group [8], which seem to be consistent with the predicted theoretical spectrum and, 
therefore, with the presence of the GZK cutoff. Meanwhile, on the other hand, we have the observations 
from the Akeno Giant Air Shower Array (AGAS A) collaboration group [9] , which reveal an abundant flux 
of incoming cosmic rays with energies above 1 x lO'^'^ eV. The appearance of these high energy events is 
opposed to the predicted GZK cutoff, and represents a great challenge that has motivated a vast amount 
of new ideas and mechanisms to explain the phenomenon [10]. If the AGASA observations are correct 
then, since there are no known active objects in our neighborhood (let us say within a radius R ~ 100 
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Mpc) able to act as sources of such energetic particles and since their arrival is mostly isotropic (without 
any privileged local source), we arc forced to conclude that these cosmic rays come from distances larger 
than 100 Mpc. This is commonly referred as the Greisen-Zatsepin-Kuzmin (GZK) anomaly. 

One of the interesting notions emerging from the possible existence of the GZK anomaly is that, since 
ultra high energy cosmic rays involve the highest energy events registered up to now. then a possible 
framework to understand and explain this phenomena could be of a quantum-gravitational nature [11, 
12, 13]. This possibility is indeed very exciting if we consider the present lack of empirical support for the 
different approaches to the problem of gravity quantization. In the context of the UHECR phenomena, all 
these different approaches motivated by different quantum gravity formulations, have usually converged 
on a common path to solve and explain the GZK anomaly: the introduction of effective models for the 
description of high energy particle propagation. These effective models, pictures of the yet unknown full 
quantum gravity theory, offer the possibility to modify conventional physics through new terms in the 
equations of motion (now effective equations of motion), leading to the eventual breakup of fundamental 
symmetries such as Lorentz invariance (expected to be preserved at the fundamental level). These Lorentz 
symmetry breaking mechanisms are usually referred as Lorentz invariance violations (LIV's) if the break 
introduce a privileged reference frame, or Lorentz invariance deformations (LID's), if such a reference 
frame is absent [14]. Its appearance on theoretical and phenomenological grounds (such as high energy 
astrophysical phenomena) has been widely studied, and offers a large and rich amount of new signatures 
that deserve attention [15]. 

To deepen the above ideas, we have adopted the loop quantum gravity (LQG) theory [16]. It is possible 
to study LQG through effective theories that take into consideration matter-gravity couplings. In this 
line, in the works of J. Alfaro et al. [17], the effects of the loop structure of space at the Planck level are 
treated scmiclassically through a coarse-grained approximation. An interesting feature of these methods 
is the explicit appearance of the Plank scale Ip and the appearance of a new length scale £ Ip (called the 
"weave" scale), such that for distances d <^ C the quantum loop structure of space is manifest, while for 
distances d > C the continuous flat geometry is regained. The presence of these two scales in thc^ effective 
theories has the consequence of introducing LIV's to the dispersion relations E = E{p) for particles with 
energy E and momentum p. It can be shown that these LIV's can significantly modify the kinematical 
conditions for a reaction to take place. For instance, as shown in detail in [18], a consequence for the 
UHECR phenomenology is that the kinematical conditions for a reaction between a primary cosmic ray 
and a CMBR photon can be modified, leading to new effects and predictions such as an abundant fiux 
of cosmic rays well beyond the GZK cutoff energy (explaining in this way the AGASA observations) . 

In this work we provide some techniques to establish and analyze new constraints on the LQG pa- 
rameters (or any LIV parameters coming from other theories). Also, we attempt to predict a modified 
UHECR spectrum consistent with the AGASA observations. The results shown here are presented in 
detail in a previous work by J. Alfaro and G. Palma [13]. This paper is organized as follows: In section 
2 we give a brief derivation of the conventional spectrum and analyze it with AGASA observations. In 
section 3 we present a short outline of loop quantum gravity and its effective description of fermionic and 
bosonic fields. In section 4 we analyze the effects of LQG corrections on the threshold conditions for the 
main reactions involved in the UHECR phenomena to take place. In section 5 we show how the modified 
kinematics can be relevant to the theoretical spectrum J{E) of cosmic rays (we will present the obtained 
modified spectrum). In section 6 we give some final remarks. 

2 Ultra High Energy Cosmic Rays 

Here we review the main steps in the derivation of the UHECR spectrum. This presentation will be useful 
and relevant to understand how LQG corrections affect the predicted flux of cosmic rays. The following 
material is mainly contained in the works of F.W. Stecker [4], Berezinsky et al. [6] and S.T. Scully & 
F.W. Stecker [7]. 

2.1 General Description 

Two simple and common assumptions used in the development of the cosmic ray spectrum are: 1) sources 
are uniformly distributed in the Universe, and 2) the generation flux F{Eg) of emitted cosmic rays from 
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the sources is well described by a power law behavior of the form F{Eg) ca Eg^" , where Eg is the energy 
of the emitted particle and 7^ is the generation index. 

One of the main quantities for the computation of the UHECR spectrum is the energy loss —E~^dE/dt. 
This quantity describes the rate at which a cosmic ray loses energy, and takes into consideration two chief 
contributions: the energy loss duo to the rcdshift attenuation and the energy loss due to collisions with 
the CMBR photons. This last contribution depends, at the same time, on the cross sections a and 
the inelasticities K of the interactions produced during the propagation of protons in the extragalactic 
medium, as well as on the CMBR spectrum. The most important reactions taking place in the description 
of proton's propagation are the pair creation and the photo-pion production: 

p + 7^p + e^+e^, and p + 7 ^ p + tt. (1) 

The photo-pion production happens through several channels (for example the baryonic A and N , and 
mesonic p and oj resonance channels, just to mention some of them) and is the main responsible in the 
appearance of the GZK cutoff. 



2.2 Some Kinematics 

To study the interaction between protons and the CMBR, it is useful to distinguish between three 
reference systems; the laboratory system K. (which we identify with the Friedman Robertson Walker 
(FRW) co-moving reference system), the center of mass (cm.) system /C*, and the system where the 
proton is at rest JC' . In terms of these systems, the photon energy will be expressed as w in /C and as e in 
K,' . The relation between both quantities is simply e = 7a;(l — (5cos9), where 7 = E/rUp is the Lorentz 
factor relating IC and IC' , E and nip are the energy and mass of the incident proton, /? = -^/l — 7~^, and 
is the angle between the momenta of the photon and the proton measured in the laboratory system 
/C. To determine the total energy E*^^ = E* + e* in the cm. system, it is enough to use the invariant 
energy squared s = E^^^. — p^^j (where Etot = E + uj and ptot are the total energy and momentum in the 
laboratory system). In this way, we have E^^ = s = nip + 2mj,e. As a consequence, the Lorentz factor 
7c which relates the /C reference system with the /C* system, is 

^ E + iv E 

Vi - (m2-f 2mpe)i/2- ^ ^ 

Let us consider the relevant case in which the reaction between the proton and the CMBR photon is 
of the type 

p + ^^a + b, (3) 

where a and b are two final particles of the collision. The final energies of these particles are easily 
determined by the conservation of energy-momentum. In the /C* system these are 

1 

Transforming this quantity to the laboratory system, and averaging with respect to the angle between 
the directions of the final momenta, it is possible to find that the final average energy of a (or b) in the 
laboratory system is 

+ (5) 

The inelasticity K of the reaction is defined as the average fractional difference K = AE/E, where 
AE = E — Ef is the difference between the initial energy E and final energy Ef of the proton (in a single 
collision with the CMBR photons) . For the particular case of the emission of an arbitrary particle a (that 

is to say p + ^ ^ p + a), expression (5) allows us to write 

(6) 

where Ka is the inelasticity of the described process. This is one of the main quantities involved in the 
study of the UHECR spectrum. 



K,b=—is + mlt,-mlJ. (4) 




3 



2.3 Mean Life t{E) 



To derive the UHECR spectrum it is imperative to know the mean hfe t{E) of the cosmic ray (or proton) 
with energy E propagating in space, due to the attenuation of its energy by the interactions with the 
CMBR photons. The mean hfe t{E) is defined through the relation t{E)~^ = {—E~^dE/dt)^^^, where 
the label "col" refers to the fact that the energy loss is due to the collisions with the CMBR photons. It 
is possible to show that the mean life t{E) can be written in the form 

= dea{e)K{e)eHl - e^/^^'^^], (7) 

where K{t) = /S.E/E is the inelasticity (with /S.E the difference between the initial and final energies of 
the proton before and after each collision), a{e) is the scattering cross section with the CMBR photons, 
KT is the temperature of the CMBR bath, 7 the lorentz factor, and eth is the threshold photon energy 
for reactions to take place. 



2.4 Energy Loss and Spectrum 

The energy loss suffered by a very energetic proton during its journey, from a distant source to our 
detectors, is not only produced by the collisions that it has with CMBR at a particular epoch. There is 
also a decrease in its energy due to the redshift attenuation produced by the expansion of the Universe. 
At the same time, such expansion will affect the collision rate through the attenuation of the photon 
gas density, which can be understood as a cooling of the CMBR through the relation T = (1 + z)Tq, 
where z is the redshift and To is the temperature of the background at the present time. To calculate 
the spectrum we need to consider the rate of energy loss during any epoch z of the Universe. 

For the present discussion, we consider a flat Friedman Robertson Walker (FRW) universe dominated 
by matter. The above assumptions give rise to the following relation between the time coordinate t and 
the redshift z: dt = —dz/Ho{l + z)^/^, where Hq is the Hubble constant at present time. Since the 
momentum of a free particle in a FRW space behaves as p a (1 + z), we will have, with the additional 
consideration m (where m is the particle mass), that the energy loss due to redshift is 

On the other hand, the energy loss due to collisions with the CMBR evolves as the background tem- 
perature changes (recall T = (1 + z)To). This evolution can be parameterized through z and is given 

by 



(9) 



Thus, the total energy loss can be expressed considering both contributions (using z instead of t): 

Equation (10) can be integrated numerically to provide the energy Eg(E. z) of a proton generated by the 
source in a z epoch and that will be detected with a energy E here on Earth. Let us express this solution 
formally by: Eg{E, z) = X{E, z)E. 

It is also possible to manipulate equation (10) to obtain an expression for the dilatation of the 
energy interval dEg/dE. To accomplish this it is necessary to integrate (10) with respect to z and then 
differentiate it with respect to E to obtain an integral equation for dEg/dE. The solution of such an 
equation is found to be 



dEgiza) , r r' dz 

^ = (l + z,)exp[y^ -(1 

where E' = {\^ z)X{E, z)E. 



+ ' dE' 



(11) 
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The total flux dj{E) of emitted particles from a volume element dV = R^{z)r'^drdfl, in the epoch z 
and coordinate r, measured from Earth at present time with energy E is 

d3{E)dE- (1 + ^)4^^2,2 . (12) 

where j{E) is the particle flux per energy, F{Eq. z)dE{) the emittcid particle flux within the range {E{). E{) + 
dEg), and n{z) the density of sources in z. As previously mentioned, it is convenient to study the emission 
flux with a power law spectrum of the type F{E) oc E~^<i . It can be shown that with such assumption, the 
relation between the emission flux and the total luminosity Lp of the source is F(E) = (7^ — 2)LpE~^s , 
To describe the evolution of the sources we shall also use a power law behavior. This will be done through 
the relation 

Lp(z)n(z) = (l + z)(-^+™)/:o, (13) 

where Cq = iyp(0)n(0). In this way, m = corresponds to the case in which sources do not evolve. If we 
consider that i?o = (1 + z)R{z) and R{z)dr = dt for flat spaces (and t; ~ 1 for very energetic particles) 
and integrating (12) from 2 = to some z = Zmax for which sources are not relevant for the phenomena, 
it is possible to obtain 

m = (7, - 2)^|^S-^' J^z^il + z,r-'/'X-^^{E, z/-^^. (14) 

The above expression constitutes the spectrum of UHECR. The volumetric luminosity Lq and the and 
m indexes are free parameters that must be fixed observationally. 



2.5 Ultra High Energy Cosmic Rays Spectrum 

To accomplish the computation of the theoretical spectrum we need information about the dynamical 
processes taking place in the propagation of protons along the CMBR. As we already emphasized, the 
most important reac:tions taking place in the description of a proton's propagation are, the pair creation 
p + 7— >p + e~+e'^, and the photo-pion production p + 7 p + -k. This last reaction is mediated by 
several channels. The main channels are N + n, A + tt, R, N + p(770), N + cj(782). The total cross 
sections and inelasticities of these processes are well known and can be used in (7) to compute the main 
time life of protons as a function of their energy. Then, with the help of expressions (11) and (14), we can 
finally find the predicted spectrum for UHECR. Fig.l shows the obtained spectrum J{E) of UHECR and 
the AGASA observed data. Again, we have selected for the theoretical spectrum J{E) the idealized case 
Emax = 00. To reconcile the data of the low energy region {E < 4 x 10^^ eV), where the pair creation 
dominates the energy loss, it is necessary to have a generation index = 2.7 (with the additional 
supposition that sources do not evolve) and a volumetric luminosity Cq = 4.7 x lO^^ergs/Mpe'^yr. It 
can be seen that for events with energies E > 4 x 10^^ eV, where the energy loss is dominated by the 
photo-pion production, the predicted spectrum does not fit the data well. To have a statistical sense 
of the discrepancy between observation and theory, we can calculate the Poisson probability P of an 
excess in the five highest energy bins. This is P = 1.1 x 10~^. Another statistical measure is provided 
by the Poisson [19]- Computing this quantity for the eight highest energy bins, we obtain = 29. 
These quantities show how far the AGASA measurements are from the theoretical prediction given by 
the curve of Fig.l. Other more sophisticated models have also been analyzed in detail [6], nevertheless, it 
has turned out that the conventional standard model of physics does not have the capacity to reproduce 
the observations from the AGASA collaboration group in a satisfactory way. 



3 Loop Quantum Gravity 

Loop quantum gravity is a canonical approach to the problem of gravity quantization. It is based on 
the construction of a spin network basis, labelled by graphs embedded in a three dimensional insertion 
S in space-time. A consequence of this approach is that the quantum structure of space-time will be of 
a polymer- like nature, highly manifested in phenomena involving the Planck scale lp. 
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Figure 1: UHECR spectrum and AGASA observations. The figure shows the UHECR spectrum J{E) 
multipUed by E^, for uniform distributed sources, without evolution, and with a maximum generation 
energy i?max = oo. Also shown are the AGASA observed events. The best fit for the low energy sector 
(JE; < 4 X 10^^ eV) corresponds to 7^ = 2.7. 

The former brief description of loop quantum gravity allows us to realize how complicated a full 
treatment of a physical phenomena could be when the quantum nature of gravity is considered, even for 
a flat geometry. It is possible, however, to introduce a loop state which approximates a flat 3-metric on 
E at length scales greater than the length scale £ ^ Ip. For pure gravity, this state is referred to as 
the weave state \W), and the length scale C as the weave scale. A flat weave \ W) will be characterized 
by C in such a way that for distances d <^ C the quantum loop structure of space is manifest, while for 
distances d > C the continuous flat geometry is regained. With this approach, for instance, the metric 
operator qat satisfies 

{W\qab\W) =5ab+0{lp/C). (15) 

A generalization of the former idea, to include matter fields, is also possible. In this case, the loop 
state represents a matter field t/j coupled to gravity. Such a state is denoted by \ W, tp) and, again, is simply 
referred to as the weave. As before, it will be characterized by the weave scale C and the hamiltonian 
operators are expected to fulfill a relation analog to (15), that is, we shall be able to define an effective 
hamiltonian such that 

H^ = {W,^\H^\W,i^). (16) 

An approach to this task has been performed by J. Alfaro et al. [17] for 1/2-spin fermions and the 
electromagnetic field. In this approach the effects of the loop structure of space at the Planck level are 
treated semiclassically through a coarse-grained approximation [21]. This method leads to the natural 
appearance of LIV's in the equations of motion derived from the effective hamiltonian. The key feature 
here is that the effective hamiltonian is constructed from expectation values of dynamical quantities from 
both the matter fields and the gravitational field. In this way, when a flat weave is considered, the 
expectation values of the gravitational part will appear in the equations of motion for the matter fields 
in the form of coefficients with dependence in both scales, C and Ip. When a flat geometry is considered, 
the expectation values can be interpreted as vacuum expectation values for the considered matter fields. 

A significant discussion is whether the Lorentz symmetry is present in the full LQG theory or not 
[22]. For the present work, we shall assume that Lorentz symmetry is indeed present in the full LQG 
theory. This assumption, jointly with the consideration that the new corrective coefficients are vacuum 
expectation values, leads us to consider that the Lorentz symmetry is spontaneously broken in the effective 
theory level. 

Specially important in the development of the effective equations of motion for matter fields is that 
they are only valid in a homogenous and isotropic system. From the point of view of a spontaneous 



6 



symmetry breakup such a system is unique and, therefore, a privileged reference frame). It is possible 
then to put the equations of motion (and therefore the dispersion relations) in a covariant form through 
the introduction of a four- velocity vector explicitly denoting the existence of a preferred system. From the 
cosmological point of view, such a privileged system does exist, and corresponds to the CMBR co- moving 
reference system. For that reason, we shall assume that the preferred system denoted by the presence of 
LIV's is the same CMBR co-moving reference frame, and will use it as the laboratory system. 

In what follows we briefly summarize the obtained dispersion relations for both, fermions and bosons. 

3.1 Modified Dispersion Relations 

The dispersion relation for fermions is found to be [17]: 

£;2 = p2 ^ 2ap2 -I- ± 2Ap -I- m^ (17) 

where the ± signs correspond to the helicity state of the described particle (note that these signs are 

produced by parity violation terms in the equations of motion). Additionally, we have defined the 
parameters a, t] and A in such a way that they depend on the scales £ and Ip in the following way: 

a = Ka{lp/CY, r] = Kr,ll, and A = lp/2£^, (18) 

where aCq,, Kr/ and k\ are adimensional parameters of order 1. 

In the case of bosonic particles the dispersion relation consists of 

= + 2ap^ + r]p^ + , (19) 

where again we have defined parameters a and r] as: 

a = Ka(/p/'C)^, and jj = Kr/ll (20) 

where Ka, Kri and k\ are adimensional parameters of order 1. 

As shown in detail in [13], only the a correction will be relevant in the computation of the spectrum. 
Therefore, in the rest of this work we shall only consider the following dispersion relation for both, 
fermions and bosons: 

E"^ =p'^ + 2ap'^ + m^. (21) 

To conclude, let us mention that the dispersion relation (21) will be used for the physical description of 
electrons, protons, neutrons, A and A'^ baryonic resonances, as well as for mesons tt, p and w. 



4 Threshold Conditions 

A useful discussion around the effects that LIV's can have on the propagation of UHECR can be raised 
through the study of the threshold conditions for the reactions to take place [18]. To simplify our 
subsequent discussions, let us use the following notation for the modified dispersion relations 

£2=p2^/(p)+m2, (22) 

where f{p) is the deformation function of the momentum p. 

A decay reaction is kinematically allowed when, for a given value of the total momentum po = 
SinitiaiP ~ SfinaiP' '^^^ ^^'^ ^ total energy value Eq such that Eq > Emin- Here E'min is the 
minimum value attainable by the total energy of the decaying products for a given total momentum po. 
To find ii^inin, it is enough to take the individual decay product momenta to be collinear with respect to 
the total momentum po and with the same direction. To see this, we can variate Eq with the appropriate 
restrictions 

Eo = ^E,{p,)+Cj{p^-^pi), (23) 

i i 
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where are Lagrange multipliers, the i index specifies the ith particle and the j index the jth vectorial 
component of the different quantities. Doing the variation, we obtain 

BE- 

-^^<=ir (24) 

That is to say, the velocities of all the final produced particles must be equal to ^. Since the dispersion 

relations that we are treating are monotonously increasing in the range of momenta p > A, the momenta 
can be taken as being coUincar and with the same direction of the initial quantity po- 

In this work, we will focus on those cases in which two particles (say a and 6) collide to subsequently 
decay in the aforementioned final states. For the present discussion, particles a and h have momenta pa 
and respectively, and the total momentum of the system is po- It is easy to see from the dispersion 
relations that we are considering, that the total energy of the system will depend only on = \pa \ and 
Vh = \Vh\- Therefore, to obtain the threshold condition for the mentioned kind of process, we must find 
the maximum possible total energy i^max of the initial configuration, given the knowledge of Po and pf,. 
To accomplish this, let us fix p^ and variate the incoming direction of pi, = npb in 

= K(P0-P6n) + E6(pb) + x(n'-l). (25) 

Varying (25) with respect to n (x is a Lagrange multiplier), we find 

= <P6/2x. (26) 

In this way we obtain two extremal situations x = ±i'ap;,/2, or simply ??.' = iu^/wa. A simple inspection 
shows that for the dispersion relations that we are considering, the maximum energy is given by = 
—v\lva, or in other words, when a frontal collision takes place. 

Summarizing the threshold condition for a two particle (a and 6) collision and subsequent decay, can 
be expressed through the following requirements: 

E^^E^>Y,E^, (27) 

final 

with all final particles having the same velocity {vi = vj for any to final particles i and j), and 

Pa-Pb = ^Pf, (28) 

final 

where the sign of the momenta J^dnaiPf given by the direction of the highest momentum magnitude 
of the initial particles. 

Our interest in the next subsections is the study of the reactions involved in high energy cosmic ray 
phenomena through the threshold conditions. To accomplish this goal through simple expressions that 
are easy-to-manipulate, we shall further use, for the equal velocities condition, the simplification 

EbrUa = EaTUb, (29) 

valid for the study of parameters coming from the region /(p) <C m^. This simplification will allow the 
achievement of bounds over the order of magnitude of the different parameters involved in the modified 
dispersion relations, which are precisely our main concern. 

In the following subsections we study the kinematical effects of LIV's through the threshold conditions 
for the reactions involved in the propagation of UHECR. Since, in this phenomena, photons are present 
in the form of low energy particles (the soft photons of the CMBR), the LQG corrections in the electro- 
magnetic sector of the theory can be ignored. LQG corrections to the electromagnetic sector, however, 
have already been studied for other high energy reactions such as the Mkn 501 7-rays [12]. Finally, let 
us recall that in the complete treatment of threshold conditions, it is possible to learn that only the a 
coefficients will be significant [13]. 
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4.1 Photo-Pion Production 7 + p ^ p + tt 

Let us begin with the photo-pion production 7 + p ^ p + tt. Considering the corrections provided in the 
dispersion relation (21) for fermions and bosons, we note that, for the photo-pion production to proceed, 
the following condition must be satisfied 

2foE^+4j;.a;> ^'('"^^ + "^-\ (30) 



where E^^ is the produced pion energy and 5a = ap — 



Too + WIt, 



4.2 Resonant Production 7 + p — > A 

The main channel involved in the photo-pion production is the resonant production of the A(1232). It 
can be shown that the threshold condition for the resonant A(1232) decay reaction to occur, is 

25aE^ +AujE>m\-ml, (31) 

where E is the incident proton energy and 6a = Up — a/\. 



4.3 Pair Creation 7 + p— >p + e+ + e 

Pair creation, 7 + p — > p + e+ + e~, is greatly abundant in the sector previous to the GZK limit. When 
the dispersion relations for fermions are considered for both protons and electrons, it is possible to find 

5a— E"^ + Eoj> me{mp + rUe), (32) 
Trip 

where E is the incident proton energy and 5a = Up — ag. 



4.4 Bounds 

In order to study the threshold conditions (30), (31) and (32), in the context of the GZK anomaly, we 
must establish some criteria. Firstly, as we have seen in section 2, the conventionally obtained theoretical 
spectrum provides a very good description of the phenomena up to an energy ~ 4 x 10^^ eV. The main 
reaction taking place in this well described region is the pair creation j + p^p + e~^ + e~ and, therefore, 
no modifications are present for this reaction up to ~ 4 x 10^^ eV. As a consequence, and since threshold 
conditions offer a measure of how modified kinematics is, we will require that the threshold condition 
(32) for pair creation not be substantially altered by the new corrective terms. 

Secondly, we have the GZK anomaly itself, which we are committed to explain. Since for energies 
greater than ^ 8 x 10^^ eV the conventional theoretical spectrum does not fit the experimental data 
well, we shall require that LQG corrections be able to offer a violation of the GZK-cutoff. The dominant 
reaction in the violated E > 8 x 10^^ region is the photo-pion production and, therefore, we further 
require that the new corrective terms present in the kinematical calculations be able to shift the threshold 
significantly to preclude the reaction. 

We begin our analysis with the threshold condition for pair production. In this case we have: 

Tfl 

5a —E^ +Euj> TTle{Tnp + Trie) , (33) 
rrip 

with 5a = Up — a^. As is clear from the above condition, the minimum soft-photon energy Wmin for the 
pair production to occur, is 

ojjnin = -^{mp + me) - 5a— E. (34) 

Hi TTlp 

It follows therefore that the condition for a significant increase or decrease in the threshold energy for pair 
production becomes \5a\ > mp{Tnp + me)/E'^. In this way, if we do not want kinematics to be modified 
up to a reference energy E^ef = 3 x 10^^, we must impose the following constraint 

lap -ae|<^^^^^^^^^= 9.8x10-22. (35) 
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Similar treatments can be found for the analysis of other astrophysical signals like the Mkn 501 7-rays 
[23], when the abscnc;c of anomalies is considered. 

Let us now consider the threshold condition for the photo-pion production. Taking only the a cor- 
rection, we have 

2foE^+4j;.a;> "^'('"^^ + "^-\ (36) 
nip + m,r 

It is possible to find that for the above condition to be violated for all energies Et^ of the emerging pion, 
and therefore no reaction to take place, the following inequality must hold 

where = KT = 2.35 x 10"** eV is the thermal CMBR energy. If we repeat these steps for the A(1232) 
resonant decay, we obtain the following condition 

OA - ap > ^ = 1.7 X 10-2^ [uj/ujof . (38) 

— 

To estimate a range for the weave scale £, let us use as a reference energy Wj-of = ujmin, where 
i-i^min is the minimum energy for the reaction to take place, in inequality (36), when the condition for 
a significant increase in the threshold condition is taken into account (for a primordial proton reference 
energy E'rof = 2 x 10^", this is LUmin ~ 2.9 x ujq), and join the results deduced from the mentioned 
requirements. Assuming that the Ka parameters are of order 1, as well as the difference between them 
for different particles, we can estimate — for the weave scale C — the preferred range 

2.6 x 10"*^ eV"* < £ < 1.6 X 10"*^ eV'S (39) 

where the lefthand and righthand sides come from bounds (35) and (37) respectively (since the A(1232) 
is just one channel of the photo-pion production, we shall not consider it to set any bound). 

If no GZK anomaly is confirmed in future experimental observations, then we should state a stronger 
bound for the difference — cxp. Using the same assumptions to set the restriction (35) when the 
primordial proton reference energy is Erei = 2 x 10^° eV, it is possible to find 

- ap\ < 2.3 x 10"^^ (40) 

In terms of the length scale C, this last bound may be read as 

£ > 1.7 x 10~" eV"\ (41) 
which is a stronger bound over £ than (35), offered by pair creation. 



5 Modified Spectrum 

In this section we show the way in which the LQG correction a affects the prediction for the theoretical 
cosmic ray spectrum. Our approach will be centered on the supposition that the LQG corrections to 
the main quantities for the calculation — such as cross sections and inelasticities of processes — are, in a 
first instance, kinematical corrections, and that the Lorentz symmetry is spontaneously broken. These 
assumptions will allow us to introduce the adequate corrections when a modified dispersion relation is 
known. 



5.1 Kinematics 

Having a spontaneous Lorentz symmetry breaking we can use the still valid Lorentz transformations 
to express physical quantities observed in one reference system, in another one. This is possible since, 
under a spontaneous symmetry breaking, the group representations of the broken group preserves its 
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transformation properties. In particular, it will be possible to relate the observed 4-niomenta in different 
reference systems through the usual rule 



= A/P^, (42) 

where = {—E,p) is an arbitrary 4-momentum expressed in a given reference system /C, is the same 
vector expressed in another given system /C', and A^" is the usual Lorentz transformation connecting 
both systems. Such a transformation will keep invariant the scalar product 

p^p^ = -E^+p^ (43) 

as well as any other product. Let us illustrate, for transformation (42), the situation in which JC' is a 
reference system with the same orientation of K. and which represents an observer with velocity /3 with 
respect to /C. In this case shall correspond to a boost in the — /3/|/3| direction, and expression (42) 
will be reduced to 

E'=j{E-p-p}, and p' = 7(p-/3E), (44) 

where 7 = (1 - /32)-i/2_ A particular case of this transformation will be that in which j3 has the same 
direction as p, and K.' corresponds to the cm. reference system, that is to say, the system in which 
p' = 0. In such a case we will have (3 = p/E and 7 = (iJ^ — p^)~^/^, jointly with the relation 

E' = Eh={E'^ -p^f/^. (45) 

In other words, the cm. energy of a particle with energy E and momentum p in JC will correspond to 
the invariant {E"^ —p^Y^"^- Furthermore, such energy is the minimum measurable energy by an arbitrary 
observer; this can be confirmed by solving equation dE' /d(3 = from the relation (44) and by verifying 
that the solution is /3 = p/E. This allows us to interpret E' = (E"^ —p^Y^"^ as the rest energy of the given 
particle. To simplify the notation and the ensuing discussions, let us introduce the variable s = {E'Y, 
where E' is given by (45). 

So far in our analysis, relativistic kinematics has not been modified. Nevertheless, a difference with 
the conventional kinematical frame is that in the present theory the product (43) will not be independent 
of the particle's energy; conversely, we will have the general expression 

p^p^, = -fa{E,p)-ml, (46) 

where fa{E,p) is a function of the energy and the momentum, that represents the LIV provided by the 
LQG effective theories. Let us note that expression (46) is just the modified dispersion relation 

E^^p^ + !a{E,p)+ml. (47) 

To be consistent fa{E,p) must be invariant under Lorentz transformations and, therefore, can be written 
as a scalar function of the energy and the momentum. 

We have already made mention of the fact that LIV's inevitably introduce the appearance of a 
privileged system; in the present discussion we will choose as such a system the isotropic system (which 
by assumption is the co-moving CMBR system), and will express fa{E,p) in terms of E and p measured 
in that system. As may be expected in this situation, fa will be a function only of the energy E and the 
momentum norm p = \p\, since no trace of a vectorial field could be allowed when isotropy is imposed. 
For example, in the particular case of the dispersion relation for a fermion, the function fa{E,p) depends 
uniquely on the momentum, and can be written as fa{p) = '^O-aV^ ■ For simplicity, we shall continue using 
/a(p) instead of fa{E,p). 

Through the recently introduced notation and the use of expression (45), the cm. energy of an a 
particle with mass nia and deformation fa{p) will be 

= Vfa{p)+ml. (48) 

Of course, the validity of this interpretation will be subordinate to those cases in which Sa = fa {p) + m-a > 
0, or, equivalent, to those states with a time-like 4-momentum. In the converse, particles with energies 
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and corrections such that Sa = fa{p) < will be described by light-like physical states if the equality 
holds, or spacc-likc physical states if the inequality holds. 

A new effect provided by LIV's is that, if a reference system where p = exists, then in that system 
the particle will not be generally at rest. To understand this it is sufficient to verify that in general the 
velocity follows v = dE/dp ^ p/E, and therefore docs not vanish at p = 0. Returning to the equations 
in (44), we can see that when (} = v = dE/dp, the following result is produced 

dE' fdE \ dp „ 

where 7^ = (1 — t;^)~^/^. Result (49) shows that the velocity of the system where the particle is at rest is 
effectively v. There emerges, then, an important distinction between the phase velocity (5 = p/E oi the 
cm. system of a particle, and the group velocity v = dE/dp of the same particle. 

The above results, for a single particle, can be easily generalized to a system of many particles. For 
instance, the total 4- momentum of a system of many particles p^'^^ = p^ will transform through the 
rule (42), and the scalar product {p^Pfi)tot will be an invariant under Lorcntz transformations (as well as 
any other product). As in the case of individual particles, we define ^/s as the total rest energy measured 
in the system of cm. That is to say: s = E^^^ —Ptot- the case in which we have a system composed of a 
proton with energy E and momentum p, and a photon (from the CMBR) with energy uj and momentum 
k (all these quantities are measured in the laboratory isotropic system /C), the s quantity shall acquire 
the form 

s = {E + ojf - (p + kf 

= {E' + uj'f-{p' + kf, (50) 

where the E' and p' quantities are measured in an arbitrary reference system. In particular we are 
interested in the system where the proton momentum p' is null; that is to say, the system in which 
E' = . If e is the photon energy in such a system, then 

s = (\/s^+e)^-e^ 

= 2y/s^ e + Sp, (51) 

where we have used the dispersion relation cj = k (or e = ') for the CMBR photons. These results will 
allow us to express the main kinematical quantities in term of e and Sp = fp{p) + rrip . For example, the 
Lorentz factor that connects the /C system with the /C' system where p' = 0, will be 

-l = E/^p. (52) 

Meanwhile, the Lorentz factor connecting K with the cm. system (that in which p' + fc = 0), will be 

£^ + w E , , 

7c = 5 7. • (53) 

Sp + + e 

As a last comment, let us note that to the first order in the expansion of the dispersion relations in 
terms of the scales £ and Ip, when we consider high energy processes such that p^ ^ f{p) + we can 
freely interchange the momentum p by the energy E in the deviation function f{p). That is to say, we 
may consider as a valid relation the following expression 

E^=p^ + f{E)+m'', (54) 

where we have made the replacement f{p) f{E). This procedure will greatly simplify the next 
discussions. 

5.2 Modified Inelsisticity: p + 7 — p + a; 

Following the same methods of section 2, let us obtain a modified inelasticity K for a process of the type 
p + 7 ^ p + a;, where x is an emitted particle that, in the present physical problem in which we are 
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interested, can be a tt, p or w meson. We note that the dispersion relation for the emerging proton (after 
a colUsion with a photon) can be written in the form: 

E'p-pl = fp{Ep)+ml, (55) 

where Ep is the final proton energy. Since the left side of (55) is invariant under Lorentz transformations, 
we can write 

iE;f-iplf = fp{Ep)+ml, (56) 

where the * denotes the quantities measured in the cm. system. On the other hand, in such a system, 
the following conservation relations of energy and momentum are satisfied: 

E;+E: = yG and (p;)^ = (p;)^. (57) 

Substituting both quantities in relation (55), wc obtain 

2^E; = s + Sp{Ep)-s^{E^). (58) 

In the same way, we also have the energy conservation relation in the laboratory system: Ep + E^ = -Etot- 
Using the definition for the inelasticity = AE/E for a process, where AE = Ei — Ef ~ E^^t — Ef, it 
is possible to rewrite the energy conservation relation in the laboratory system in terms of through 
expressions E.^ = K^E, and Ep = {1 — Kx)E, where E is the initial energy of the initial proton. Having 
done this, equation (58) now acquires the form: 

2V^E; = s + Sp[{l - K^)E] - s^[K^E]. (59) 

To simplify the development of the inelasticity, let us write the former relation as E* = F{E, Kx), where 
F = F{E, Kx) is defined through 

" 1 {s + Sp[{l-Kx)E]-Sx[KxE]). (60) 



On the other side, the Lorentz transformation rules give us the relation between the proton energies in the 



laboratory system and the cm. system respectively. This relation is -Ep/ 7c = Ep+Pcy E*"^ — Sp(Ep) cosO. 
Joining this expression with (60), it is possible to find the general equation for K^,: 

(1 - Kx)V^ = (^F{E, Kx) + ^F^E, Kx) - Sp[{l - Kx)E] cosO) . (61) 

It should be noted, however, that the solution for K^ from (61) will depend in the 6 angle. For this 
reason, once this last equation is resolved, it is convenient to define the total inelasticity K as the average 
of Kx with respect to the 6 angle. That is to say 



K=- I Kxde. (62) 

It is relevant to mention that now, as opposed to result (6), the inelasticity K will be a function of both, 
the energy E of the initial proton and the energy e of the CMBR photon. 



5.3 The — ^ Sa = + fa{E) Prescription 

Let us recall our interpretation relative to the fact that = {fa{Ea) + i^a)^^^ can be understood as 
the rest energy of a particle a, as a function of the energy Ea that it has in the laboratory system /C. As 
we have already emphasized, such an interpretation will be valid for particles with time like 4-momenta. 

In the reactions given between high energy protons and the photons of the CMBR, the whole scenario 
consists of the collision between two particles, p and 7, with the subsequent production of a certain 
number of final particles. Let us suppose that a is one of these particles in the final state. The knowledge 



13 



of the inelasticity K for the reaction will allow us to estimate the average energy {Eg) with which such a 
particle emerges (since K provides the average fraction of energy with which such a particle is produced). 
That is to say, on average, the rest energy of the final particle a will be s]J'^ = [fa{{Ea)) + m^]^/^. 
Moreover, the knowledge of the inelasticity K will allow us to express Sa as a function of the energy E 
of the initial proton: 

Sa = Sa{E). (63) 

Following our previous interpretation, we can view the recently described process as a reaction between a 

1 /2 . . . . . 1 /2 

proton with mass Sp , which loses energy emitting particles a with mass Sa calculated in the previous 
form. This idealized reasoning gives us a clear prescription to kinematically modify those dynamical 
quantities with which we must work and where energy conservation is involved. This prescription is: 

ml^Sa{E) = UE)+ml, (64) 

where we have expressed correction fa as a function of the initial energy of the incident proton. 

Prescription (64) establishes the notion of an effective mass which is dependent on the initial energetic 
content of a reaction. As a consequence, given the explicit knowledge of the dependence that a cross 
section has on the masses and energies of the involved states, to obtain the modified version, it will be 
appropriate to use the discussed prescription. 



5.4 Redshift 

Another important problem related to the introduction of LIV's in the dispersion relations is whether the 
redshift relation for the propagation of particles in a FRW universe is modified. This could be of great 
relevance because of the large distances involved in cosmic ray propagation and, therefore, the possible 
cumulative effects. In the case of a FRW universe with scale factor R, it is possible to find that the 
propagation of particles with momentum = g'^^PiPj (with gij the spatial part of the FRW metric) 
obeys the following equation 

^{PR) = 0, (65) 

where r is a proper parameter describing the path followed by the particle. This result is regardless of 
any particular dispersion relation. 



5.5 Spectrum and Results 

Introducing the above modifications to the different quantities involved in the propagation of protons (like 
the cross section a and inelasticity K), we are able to find a modified version for the UHECR energy loss 
due to collisions. Since the only relevant correction for the GZK anomaly is a, we focused our analysis 
on the particular case f{p) = 2ap'^. To simplify our model we restricted our treatment to the case a > 
(consistent with the effective mass interpretation) and used only am ^ 0, where Um is assumed to have 
the same value for mesons tt, p and cj. 

Fig. 2 shows the modified energy loss t{E) for UHECR obtained for different values of am- It can 
be seen therefore how the corrections can affect the main life time of protons propagating through the 
CMBR, allowing a strong improvement in the distances that protons can reach before loosing their 
characteristic energy (for energies greater than 1 x 10^*^ eV). The effects that the LQG corrections have 
on the propagation of UHECR are manifest through a decay of the energy loss in the range _B ^ 1 x 10^*^ 
eV. To understand this, recall relation (36) for the threshold condition of photo-pion production: 

2foE^+4j;.c.> "^'('"^^ + "'-) . (66) 
nip + rriTT 

As we saw in section 4, the condition for a significant increase or decrease in the energy threshold can be 
calculated as \Sa\ > {2mp + mjr)(mp + m^)/2£J^. Therefore, for a given value of 6a > 0, the energy at 
which the LIV effects start to take place is 

-B^ = + m^){mp + m„). (67) 
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10" 10'» lO^o 10^1 

Energy E [eV] 

Figure 2: Modified energy loss for UHECR due to collisions. The figure shows the case a™ ^ 0, for 
three different values of the weave scale L. Curve 1: a™ = 9 x 10^^'^ (£ ~ 8.6 x 10^^^ eV~^); curve 2: 
a„i 5 X 10~^^ (£ ~ 1.2 X 10"^'' eV~^); curve 3: a™ ~ (without modifications). 

In the case am = 9 x 10"^'^ (curve 1 of Fig. 2), this energy \s, E — 1.1 x 10^" eV, while in the case 
ckm = 5 X 10~^^ (curve 2) this corresponds to i? = 1.5 x 10^" eV. Beyond these energy scales, at about 
E ^ 2 X 10^" eV, a sharp decay is observed in the behavior of the curve. This is due to the fact that the 
modified inelasticity K will strongly constraint the energy-momentum phase space accessible to the final 
states depending on the initial energy E that the primary proton carries (recall that now X is a function 
of the energy E of the incident proton and the energy e of the CMBR photon). 

Also, we can find the modified version of the UHECR spectrum for am ^ 0. Fig. 3 shows the AGASA 
observations and the predicted UHECR spectrum. The Poisson probabilities of an excess in the five 
highest energy bins for the three curves are Pi = 3.6 x 10^'*, P2 — 2.6 x lO""* and P3 — 2.3 x 10"'*. 
The Poisson for the eight highest energy bins are Xi = 10? X2 — 10-9 and X3 = 11-2 respectively. 
The possibility of reconciling the data with finite maximum generation energies is significant given that 
conventional models require infinite maximum generation energies -Emax for the best fit. For the lower 
part of the spectrum (under E = A x 10*^ eV), the parameters under consideration leave the spectrum 
completely unaffected. This is due to the fact that in such a region the dominant reaction is the pair 
production, which has not being modified to obtain the spectrum. A more accurate study on this issue 
would require the computation of a modified inelasticity for the pair creation. Meanwhile, we must 
content ourself with the scmiqualitative criteria given in section 4 to rule out the parameters. 

6 Conclusions 

We have seen how the kincmatical analysis of the different reaction taking place in the propagation of 
ultra high energy protons can set strong bounds on the parameters to the theory. In comparison with our 
previous work, wc have eliminated some previously open possibilities by the particular study of the pair 
creation j3 + 7^p + e+ + e^, in the energy region where this reaction dominates the proton's interactions 
with the CMBR. In this way, the only possibility still open is the correction a. If this is the case, a 
favored region for the scale length C estimated through the threshold analysis would be 

2.6 X 10"^* eV"* < £ < 1.6 X 10"" eV"^ 

Similarly, the kinematical corrections can be studied in more detail when their effects are considered in the 
theoretical spectrum. In this regard, we have seen how to develop a modified version of the inelasticity for 
the photo-pion production, and its implications in the mean life time of a high energy proton as well as on 
the spectrum. To accomplish this last task we have only assumed a spontaneous Lorentz symmetry break 
up in the effective equations of motion, allowing the use of Lorentz transformations on the dispersion 
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Figure 3: Modified UHECR spectrum and AGASA observations. The figure sliows tlie modified spectrum 
J{E) multiplied by E^, for uniform distributed sources and without evolution, for the case am = 1.5 x 
10~^^ (£ ~ 6.7 X 10""'^^ eV~^). Three different maximum generation energies -Emax are shown. These 
are, curve 1: 5 x 10^° eV; curve 2: 1 x 10^^ eV; and curve 3: 3 x 10^^ eV. 



relations. Therefore, result (61) can be used in a more general context than the special case offered by 
the LQG framework. 

Future experimental developments like the Auger array, the Extreme Universe Space Observatory 
(EUSO) and Orbiting Wide- Angle Light Collectors (OWL) satellite detectors, will increase the precision 
and phenomenological description of UHECR. On the more theoretical side, progress in the direction 
of a full effective theory, with a systematic method to compute any correction with a known value for 
each coefficient, is one of the next steps in the "loop" quantization programme [25]. Therefore, it is 
important to trace a phenomenological understanding of the possible effects that could arise as well as 
the constraints on LQG, in the high and low energy regimens (for other phenomenological studies of LQG 
effects, see for example [26]). 

Recently[27] the effective theory approach have been under criticism. 

In [28] , the type of terms we used to explain AGASA data are computed within the Standard model, 
depending on one arbitrary parameter. The criticism of [27] does not apply to this formalism. 
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In the following we are going to consider... 



Figure 1: Scatter plots of the combined CTD+VXD tracks: the left figure for tracks with ptiab > 
lOGeV/c, the right figure for tracks from the region (18) with ptiab > 100 GeV/c. 
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